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One of the possible low-energy consequences of string theory is the addition of a Chern-Simons 
term to the standard Einstein-Hilbert action of general relativity. It can be argued that the 
quintessence field should couple to this Chern-Simons term, and if so, it drives in the linearized 
theory a parity-violating interaction between the gravito-electric and gravitomagnetic fields. In this 
paper, the linearized spacetime for Chern-Simons gravity around a massive spinning body is found 
to include new modifications to the gravitomagnetic field that have not appeared in previous work. 
The orbits of test bodies and the precession of gyroscopes in this spacetime are calculated, leading 
to new constraints on the Chern-Simons parameter space due to current satellite experiments. 
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I. INTRODUCTION 

The study of modifications of the theory of general rel- 
ativity has been of interest ever since Einstein first formu- 
lated general relativity in 1915. Particularly interesting 
are modifications that introduce terms to the Einstein- 
Hilbert action that are second order in the curvature, as 
such modifications represent high-energy corrections to 
the Einstein-Hilbert action that might arise in quantum 
gravity. Chern-Simons gravity is an example of such a 
second-order modification of the Einstein-Hilbert action. 

Chern-Simons modifications to gravity were first con- 
sidered in 2+1 dimensions Refs. @, HJ investigated 
the structure of these theories in 3+1 dimensions and 
showed how they could arise as a low-energy consequence 
of string theory. Ref. [J] considered some early-universe 
implications of such theories. Refs. [|| investigated how 
Chern-Simons terms might participate in leptogenesis. 
Ref. [6| renewed the investigation of Chern-Simons grav- 
ity, working out the linearized equations of the theory 
and their implications for gravitational waves. Most re- 
cently, Refs. 0] solved the linearized Chern-Simons field 
equations around a collection of spinning point masses. 
In much of the work on Chern-Simons gravity, the Chern- 
Simons term is coupled to a scalar field (as detailed be- 
low), and this scalar field is assumed to be spatially ho- 
mogeneous but time varying. This assumption can be 
motivated by arguments analogous to those that have 
been made suggesting that the quintessence field should 
be coupled to the Chern-Simons term of electromag- 
netism [8|. 

Chern-Simons gravity has thus far eluded constraints 
from Solar System tests of weak-field gravity because 
it is indistinguishable from general relativity for all 
spacetimes that possess a maximally symmetric two- 
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dimensional subspace and for all conformally flat space- 
times |2j . Therefore, the Schwarzschild spacetime as well 
as the Robertson- Walker spacetime are also solutions of 
the Chern-Simons gravitational field equations. Distin- 
guishing Chern-Simons gravity from general relativity re- 
quires considerations of spacetimes that are not spheri- 
cally symmetric, such as the spacetime around a spin- 
ning body. To this end, Refs. Q investigated the Chern- 
Simons modifications to the motion of bodies around a 
spinning point mass and found that the motion was in- 
distinguishable from that in general relativity. 

In this paper we take further steps to link Chern- 
Simons gravity to current and forthcoming experimental 
tests of weak-field gravity. We assume, as in other recent 
work, that the scalar field coupled to the Chern-Simons 
term is time varying but spatially homogeneous. We then 
determine the spacetime around an extended spinning 
mass and find that it differs from the spacetime around 
a spinning point mass. We determine the orbits of test 
particles and the precession of gyroscopes moving in this 
spacetime and find that the Chern-Simons modification 
does lead to observable deviations from the predictions of 
general relativity. These deviations allow us to evaluate 
constraints to the Chern-Simons parameter space from 
current satellite experiments, as well as those regions to 
be probed with forthcoming experiments. 

The paper is organized as follows: Section II defines 
the theory and derives the gravitational field equations. 
Section III considers the linear theory and derives the 
gravitomagnetic equations of motion (the Chern-Simons 
Ampere's law). Section IV discusses the solution for the 
gravitomagnetic field around a spinning massive body. 
In Section V, we consider the orbital precession of test 
bodies in this spacetime, as well as the orbital preces- 
sion of gyroscopes, and we determine the regions of the 
Chern-Simons-gravity parameter space that are probed 
with the LAGEOS and Gravity Probe B satellites. We 
conclude briefly in Section VI. Appendix A shows how the 
Chern-Simons Lagrangian we work with may be derived 
from a string-theory action, and Appendix B outlines the 
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derivation of the gravitomagnetic vector potential around 
a spinning sphere. 



II. CHERN-SIMONS GRAVITY 



We consider the theory defined by the action 
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where £ ma t is the Lagrangian density for matter, g = 
det g^ v is the determinant of the metric g^ v , R is the 
Ricci scalar (with the convention R x fll , K = T x K + ■ ■ ■ 

for the Riemann tensor), and RR is a contraction of the 
Riemann tensor and its dual: 

RR = R^ s R a ^ s , (2) 

where the dual of the Riemann tensor is defined by 
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where e aTa p is the Levi-Civita tensor, including a factor 
of y/—g. Finally, £ is a new length scale, a parameter 
of the theory, and k 2 = 8irG, where G is Newton's con- 
stant. Throughout this paper we take Greek indices to 
range from to 3. Appendix A shows how such an action 
may arise in string theory. This action is different from 
the action considered in Ref. [6( in that here 9 is a dy- 
namical scalar field with a canonical kinetic term, so the 
£ parameter is required to make the action dimcnsionlcss. 
The equation of motion for 9 is given by 

ue =^-h mk - (4) 

The gravitational field equations take the form 
2 

Gfiv G^ v = k T^ u , (5) 



where is the Einstein tensor, T^„ is the stress-energy 
tensor for the scalar field and the matter Lagrangian, and 
we refer to as the Cotton- York tensor 1 , 
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(d a 0) (e^ IQ/3 V a R v + e™ al3 V a R^} (6) 
V T (d a 8) {R T ^ v + R Tva ^ 



Appendix A provides an alternative expression for the 
Cotton- York tensor. 

Ref. @ notes that if 9 is a non-dynamical field (a La- 
grange multiplier), the theory cannot accommodate a 
spacetime with a nonzero RR because the Cotton- York 
tensor would have a non-zero divergence. However, if 9 
is a dynamical field, then the theory can indeed accom- 
modate spacetimes with nonzero RR since we have 
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where T^ v is the stress-energy tensor for 9. We see that 
whereas the scalar-field stress-energy and the Cotton- 
York tensors are separately conserved when RR = 0, 
the divergence of the scalar field stress-energy tensor is 
precisely balanced by the divergence of the Cotton- York 
tensor for non-zero RR due to the novel coupling be- 
tween the scalar field and gravity. 



III. THE CHERN-SIMONS 
GRAVITOMAGNETIC EQUATIONS 

We begin with a perturbation to the flat metric [using 
signature (—+++)], 



(8) 



and compute the linearized Einstein and Cotton- York 
tensors, 
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1 We note that this definition differs from the usuai expression 
for the four-dimensional Cotton- York tensor (see Ref. [6|]). 



where □ is the flat-space dAlembertian and the comma 
denotes partial differentiation. Since we will require be- 
low only the gravitomagnetic fields, we will be primarily 
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interested in the time-space components of the linearized 
field equations. 

In this paper, we suppose that the scalar field depends 
only on cosmic time, 9 — 9(t), the assumption being 
that 9 is either a quintessence field or some other field 
that somehow echoes the arrow of time associated with 
the cosmic expansion. This choice implies that the field 
equations are not Lorentz invariant in the Solar System 
since d a 9 points in the cosmic time direction and couples 
to local gravity through the Cotton- York tensor [Eq. (6)]. 
We note that a nonzero RR will source spatial variations 
in 9 through Eq. (4). By restricting 9 to be spatially ho- 
mogenous, we are effectively treating 9 as a nondynamical 
field, and we leave a full dynamical treatment to future 
work. Finally, we neglect corrections due to the motion 
of the Earth with respect to the rest frame of the cosmic 
microwave background. 

We work with the trace-reversed metric perturbation, 

hpv = h^v - -rj^uh, (11) 

and impose the Lorenz-gauge condition, d^h^ v = 0, to 
obtain the linearized time-space field equations, 

o 

G linear ti ,„2 /^linear ,.2ti /i o\ 

ot ~ 3 £k c oi = ~ K T oi, ( 12 ) 

with 

G|T ar = \oh 0h (13) 

CcT ar = e -e% k d^ah k 0l (14) 

where the dot denotes differentiation with respect to time 
and Latin indices are purely spatial and range from 1 to 
3. The stress-energy tensor for 9{t) is diagonal, so it does 
not contribute to the time-space field equations. 

Let t a be a unit vector in the coordinate time direction, 
and then define the 4-vector potential of this linearized 
theory, 

A-n = --h^t" = -^V>- (15) 

We consider a source with mass density p, mass current 
J and negligible pressure, so we can express the matter 
stress-energy tensor as 

Tpy = 2t( Jl J v ) - pt^t v , (16) 

where J M = —T^t u = (— p, J). In general relativity, the 
time-space components of the linearized field equations 
take the form 

d^Ai = -4^GJ, ; , (17) 

which is (nearly) identical to Maxwell's equations for the 
vector potential in Lorenz gauge (d^A^ — 0). Given our 
definition of A 11 , the Lorenz-gauge condition for A^ is 
implied by our earlier gauge choice for h^. 



The classically 'physical' fields (i.e., those that enter 
into the geodesic equation) E and B are given by 

E* = d t Ao~d A t , (18) 
B i = e aijk djAki ( 1Q ) 

where we have defined e 0l: > k = 1. Two of the Maxwell 
equations, 

V-B = 0, (20) 

Vx£ = -_, (21) 

are a direct consequence of the way in which the E and 
B fields are defined in terms of the vector potential, and 
so these two equations will be the same in Chern-Simons 
gravity. Gauss' law, which follows from the time-time 
component of the field equation, is now 

V • E = AttG{ P + pg) (22) 

where pg is the energy density of the scalar field 6(t) and 
is uniform throughout the Solar System. Since pg cannot 
be larger than the mean cosmological energy density, it 
must be negligible compared to the density of the source 
p, and we do not consider it further. The only significant 
modification will be to Ampere's law, which, for Chern- 
Simons gravity, is now given by 

dE 1 

VxB-- OB = AttGJ, (23) 

ot m cs 

where we have defined m cs = —3/(£k 2 9). 

Given the metric perturbation represented by the grav- 
itomagnetic potential and neglecting the time variation 
of the metric, slowly moving particles travel on geodesies 
such that a 'Lorentz force law' of the form, 

a= -E -4vx B, (24) 

is obtained. Therefore, as in electrodynamics, only the 
physical fields, and not the potentials, have physical rel- 
evance. 

We furthermore note that RR can be expressed in 
terms of gravito-electric and gravitomagnetic fields as 

RR = -l&idiE^dkB^rf 1 + rf l rf k ). (25) 

Unlike the case with Maxwell fields @ , it is not sufficient 
for the fields to have a non-vanishing E ■ B in order to 
have a non-trivial coupling between gravity and the scalar 
field. The best example of a gravitational source which 
produces a non-vanishing RR is a spinning, spherical 
body. 

IV. GRAVITOMAGNETISM DUE TO A 
SPINNING SPHERE IN CHERN-SIMONS 
GRAVITY 

We are now in a position to calculate the gravitomag- 
netic field in Chern-Simons gravity for a spinning body. 
Appendix B provides details of the calculation. 
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We consider a homogeneous rotating sphere, and so 
the mass current is 



J = p[ujxr\e(R-r), 



(26) 



where R is the radius of the rotating body, p is its den- 
sity, Co is its angular velocity, r is the distance from the 
origin, and is the Heaviside step function. As detailed 
in Appendix B, the field equation, Eq. (|23| . is rewrit- 
ten as an equation for A and is solved by imposing the 
condition that the metric be continuous everywhere and 
that the gravitomagnetic field be finite and well-behaved 
at the origin; the resulting vector potential is given in 
Appendix B. We note that in deriving this solution we 
have assumed that the time derivative of m cs is negligi- 
ble. The gravitomagnetic field is then obtained by taking 
the curl of A and may be written as B = -Bgr + -Bcs> 
where 
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(27) 



j£ [2w + 3f x (f x CS)] 



x (fx to), r < R, 



r>R, 



is the gravitomagnetic field inside and outside a spinning 
sphere in general relativity, and 

B cs = 4nGpR 2 {D 1 (r)Lo + D 2 (r)fxio 

+D 3 {r)fx (fxco)}, (28) 

is the new contribution in Chern-Simons gravity. Inside 
the sphere (r < R), 



(m cs R) 2 
m cs r 



+ — y2{m cs R)j 1 (m cs r), 
r 

+ m cs Ry 2 (m cs R)ji{m cs r), 



Di(r) = 

D2ir) = (m cs R) 2 
D 3 (r) = m cs Ry 2 (m cs R)j 2 (m cs r), (29) 

and outside the sphere (r > i?) 

Di{r) = — j 2 {m cs R)yi(m cs r), 
r 

D 2 {r) = m cs Rj 2 (m cs R)yi(m cs r), 

D 3 (r) = m cs Rj 2 (m cs R)y 2 (m cs r), (30) 

where je(x) and yt{x) are spherical Bessel functions of 
the first and second kind. We see that the Chern-Simons 
terms alter the components of the gravitomagnetic field 
along the rotation axis Co and f x (f x Co), and they 
also introduce a new component perpendicular to the 
plane defined by Co and f. In other words, while in gen- 
eral relativity a toroidal mass current implies a poloidal 
gravitomagnetic field, the parity violation introduced in 
Chern-Simons gravity introduces a toroidal component 
to the gravitomagnetic field. Something similar occurs 
in Chern-Simons electromagnetism [9j, although the de- 
tailed fields differ since the W 2 B term in Eq. (|23|) is sim- 
ply B in the electromagnetic theory. 



The Chern-Simons addition to Ampere's law, Eq. (|23[) . 
changes that equation from a first-order differential equa- 
tion for B to a second-order differential equation. As a 
result, the Chern-Simons modification to the gravitomag- 
netic field cannot, in general, be obtained by perturbing 
around the general-relativistic result, as the solution in 
Eq. (|28p shows. In Chern-Simons gravity, the gravito- 
magnetic field oscillates with distance outside the source, 
and the amplitude of the oscillating field is not necessar- 
ily smaller than the general-relativistic gravitomagnetic 
field. Still, we expect from Eq. (|23[) that as m cs — > 00, the 
general-relativistic solution should be recovered. This oc- 
curs since the oscillatory terms vanish as m cs — ► 00, and 
so the effects on geodesies of these new terms will vanish. 

If we take Co to lie in the z direction, then the Chcrn- 
Simons gravitomagnetic field has a nonzero azimuthal 
component B^. Since B^ ^ 0, one cannot find a coordi- 
nate transformation that causes both A r and Ag to van- 
ish. This is at odds with claims (see, e.g., Ref. [10() that a 
metric for stationary axisymmetric spacetimes in Chern- 
Simons gravity can always be found with hte — h tr = 0. 
In general relativity, one can always find a coordinate sys- 
tem for which A r = Ag = for a stationary axisymmetric 
spacetime sourced by rotating perfect fluid. However, the 
proof of this statement assumes time-reversal invariance 
of the fundamental equations. This invariance implies 
that the metric components possess the same symmetries 
as the source, namely invariance under a transformation 
that takes t — » — t and <fi — > —<fi. In that case, A r and 
Ag must be zero to keep the line element invariant under 
the same transformation. In Chcrn-Simons gravity, time- 
reversal invariance is explicitly broken by the rolling of 
the scalar field, 9^0, and it is straightforward to ver- 
ify that our solution for A, given in Appendix B, implies 
that A r and Ag are both odd under time reversal. Con- 
sequently, the line element has the same symmetry as the 
source even though A r and Ag are nonzero. 

Inspection of our solution for the vector potential given 
in Appendix B shows that it differs from the solution 
for a point-like mass-current dipole (i.e., a gravitomag- 
netic dipole) obtained by Alexander and Yunes (AY) [3] • 
When applied to a single spinning source, the metric 
given by Refs. Q corresponds to a vector potential 
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(31) 

This vector potential is an exact solution to Eq. (f23|) out- 
side of a spinning sphere, and we can see that every term 
in Aay also appears in our solution for A. The additional 
oscillatory terms in our solution constitute a homoge- 
neous solution to Eq. (|23|) . but without these terms, A 
would not be continuous across the surface of the sphere. 
Furthermore, only these oscillating terms contribute to 
.Bcs because V x Aay = V x Aqr- The inclusion of oscil- 
latory terms results in a Chern-Simons gravitomagnetic 
field that differs from general relativity, so we may use 
observations of the motion of test bodies in the Earth's 
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gravitomagnetic field to constrain Chern-Simons gravity. 



V. ORBITAL AND GYROSCOPIC PRECESSION 

A. Orbital precession 

In order to investigate how the Chern-Simons gravit- 
omagnetic field will affect the motion of test particles 
around the Earth, we will use what are known as the 
Gaussian perturbation equations [ill, E3 • Details of how 
these equations are applied to gravitomagnetic forces are 
discussed in Ref . [l3j ; here we give only a brief introduc- 
tion. The Gaussian perturbation equations give the time 
variation of the Keplerian orbital elements in the pres- 
ence of a perturbing force. In our case we take the gravit- 
omagnetic force, — iv xB, as a small perturbing force and 
approximately solve the equations given in Ref. [HI . We 
will concentrate on analyzing the secular (non-periodic) 
time variation of the longitude of the ascending node 2 , 
f2, but note that other Keplerian elements will also vary 
due to the terms introduced by Chern-Simons gravity. 
The time variation of fl has been well studied since, in 
general relativity, it is connected with the Lense-Thimng 
drag |H, 



^GR = 



2GL 



a 3(!_ e 2)3/2> 



(32) 



where L is the magnitude of the angular momentum of 
the central body, a is the semi-major axis of the orbit of 
the test body, and e is the orbit's eccentricity. Finally, in 
order to evaluate the secular perturbations, we approxi- 
mate the orbit of the test body as circular (i.e., e = 0, a 
good approximation for current measurements), and we 
average the perturbing force over one orbital period to 
obtain 



^cs 

a 



GR 



15 -jph( m csR)yi(m c 



(33) 



where tics is the precession due to i?cs ■ The total preces- 
sion is f^GR. + ^cs- We note that tics is an even function 
of m cs . 

Recent measurements of laser ranging data to the LA- 
GEOS I and LAGEOS II satellites have measured tl to 
within 10% of its value in general relativity 16]. Re- 
quiring that the Chern-Simons contribution does not ex- 
ceed 10% of the general relativity result, we find that 
we can place a lower limit to the Chern-Simons mass, 
|m cs | > 0.001 km - , as shown in Fig. [T] 



2 The longitude of the ascending node is denned to be the angle 
between a stationary reference line and the line connecting the 
origin of the coordinate system and the point where the orbiting 
body intersects the XY reference plane as it is moving upwards 
(see Ref. Q). 



0.4 



0.2 - 



■a 



-0.2 



-0.4 - 




0.001 



0.01 



(km- 



FIG. 1: The ratio fics/ficR for the LAGEOS satellites orbit- 
ing with a semimajor axis of a ~ 12, 000 km. A 10% verifi- 
cation of general relativity [l(| (the shaded region) leads to a 
lower limit on the Chern-Simons mass of |m cs | > 0.001 km -1 . 
A 1% verification of the Lense-Thirring drag will improve this 
bound on m ca by a factor of roughly five. 



The Laser Relativity Satellite (LARES) mission pj} 
proposes to deploy a new laser ranging satellite and is 
predicted to measure £1 to within 1% of its value in gen- 
eral relativity. With this improvement the bound on m cs 
is increased by a factor of roughly five. 



B. Gyroscopic precession 

The Earth's gravitomagnetic field will also cause a pre- 
cession of gyroscopes moving in the spacetime. A gyro- 
scope will undergo precession due to two torques. One 
is known as the geodetic precession and is independent 
of the Earth's gravitomagnetic field. The other torque is 
due to a coupling to the gravitomagnetic field and results 
in a rate of change of the spin of a gyroscope given by 

M 



S = 2BxS, 



(34) 



where S is the angular momentum of the gyroscope. 

NASA's Gravity Probe B (GPB) mission is curre ntly 
attempting to measure this gyroscopic precession [l9j |. 
GPB consists of a satellite, in a polar orbit at an altitude 
of about 640 km, that contains four drag-free gyroscopes 
and a telescope. The gyroscopes are initially oriented 
such that their spins are aligned parallel to the optical 
axis of the telescope, which is pointing within the plane 
of the orbit. The telescope points towards a guide star, 
allowing a measurement of the precession of the direc- 
tion of the spins of the gyroscopes. Geodetic precession 
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results in an annual precession in the North-South di- 
rection of about 6600 milliarcseconds (mas) whereas the 
general relativistic gravitomagnetic field causes an annual 
East- West precession of around 42 mas 




0.01 0.1 
\m cs \ (km-') 

FIG. 2: The ratio < tcs/4 > GR for Gravity Probe B in a polar 
orbit at an altitude of approximately 640 km. A 10% verifica- 
tion of general relativity (the shaded region) leads to a lower 
limit on the Chern-Simons mass of |m cs | > 0.01 km -1 , an 
order of magnitude improvement over the LAGEOS result. 

With the Chern-Simons expression for the gravitomag- 
netic field, given in Eq. (f2"5)) . it is straightforward to cal- 
culate the resulting gyroscopic precession for a polar or- 
bit (applicable to GPB). Relative to the general relativity 
result, we find 

i = 15—j 2 (m cs R) [yi(m cs a) + m cs ay (m cs a)} ,(35) 

$ GR R 2 

where <f> = = 4>gr + ^cs is the rate at which 

the angle of axis <E> changes in time due to the gravito- 
magnetic field. We note that <I>cs is an even function of 
m cs . 

It was initially projected that GPB would achieve a 
percent-level measurement of the gravitomagnetic con- 
tribution to ^gr- However, since its launch in 2004, it 
has encountered several unexpected complications that 
will degrade the precision of the tests of gravity [2(| , al- 
though the extent of the degradation has yet to be re- 
ported. In Fig. [21 we plot Eq. ([55]) for a GPB detection 
of the gravitomagnetic precession to within 10% of its 
value in general relativity. 

We have idealized the Earth to be a sphere of con- 
stant density throughout this work, when in reality, it is 
an oblate spheriod with layers that have different mean 
densities. However, we expect that the non-spherical 
corrections would affect both the general relativity and 
Chern-Simons calculations similarly and, to the accuracy 



we require, are negligible when we consider the ratio be- 
tween general relativity and Chern-Simons results. Fur- 
thermore, it is easy to generalize our results to spheres 
with layered density profiles because B depends linearly 
on p. We replaced our model of a homogeneous Earth 
with a model of the core and mantle and we found that 
the amplitudes of the oscillations in f^cs and <I>cs were 
not affected. We conclude that our constraints on m cs 
are not sensitive to the details of the density profile of 
the Earth. 



VI. CONCLUSIONS 

The addition of a Chern-Simons term to the action for 
gravity is of interest as it may arise as a low-energy limit 
of string theory. The theory and formalism of this mod- 
ification of gravity have been worked out in a number of 
previous papers, and some of the early-Universe conse- 
quences of such a term have been investigated. However, 
there has been little work on tests of such modifications 
in the present Universe. 

In this paper, we have calculated the linear-theory 
spacetime around a spinning massive body, finding new 
corrections that were overlooked in previous work. The 
gravitomagnetic field in Chern-Simons gravity differs 
from that in general relativity in two ways: (1) there 
is an oscillating component, and (2) there is a toroidal 
component to the gravitomagnetic field that arises as a 
consequence of the parity-breaking nature of the theory 
and that has no counterpart in ordinary general relativ- 
ity. 

We then determined the precession of orbits of test 
particles in this spacetime and also of gyroscopes moving 
in this spacetime. We showed that current constraints 
from the LAGEOS satellites restrict the inverse Chern- 
Simons mass parameter m" 1 to be less than roughly 1000 
km, corresponding to a mass constraint m cs >2x 10~ 22 
GeV. This bound may be improved by a factor of 5-10 
by future observations. 

The mass parameter m cs is related to the more fun- 
damental parameters £ and 9 of the theory through 
m cs = —3/(8ttG£9), where £ is a length parameter that 
enters into the Chern-Simons Lagrangian, and 9 is pre- 
sumably related to the time variation of the quintessence 
field. In principle, a precise constraint to £ can be de- 
rived once the precise nature of the field (a quintessence 
field?) 9 and its time evolution are specified. We leave 
such model building for future work. 
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APPENDIX A: A STRING INSPIRED 
DERIVATION OF THE CHERN-SIMONS FIELD 
EQUATIONS 

The effective 4-D string action for heterotic and type 
II string theory can be written as @, Hl[ 



S = / d 4 :rV^5 



2k 



;R - OiHuvxH 11 



, (Al) 



where R is the Ricci scalar, H^ v \ is the Kalb-Ramond 
(KR) three-form field strength, and a is a constant with 
units of length squared. We are neglecting numerous 
terms, including Gauss-Bonnet terms, dilaton terms, and 
matter terms, some of which depend on compactification. 
The Kalb-Ramond field is written in differential-form no- 
tation as 



(A2) 



where B is a two-form field (known as the KR field) 
and ui is the Lorentz-Chern-Simons term. The Lorentz- 
Chcrn-Simons three-form can be written in terms of the 



spin connection uj as [22 1 



(A3) 



where the trace is over the suppressed vector indices of 
the spin connections. We then have the identity, 



dH = -Tr(RAR), 



(A4) 



associated with the KR field strength, where R is the 
Riemann tensor and the trace is over the tensor indices; 
the right-hand side is also known as the Hirzebruch den- 
sity. Taking the Hodge dual of the Hirzebruch density, 
we obtain 



^{RhR) = ^'' x R a0IJiV R a \ x 
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RR. (A5) 



Let us now consider the equation of motion for the 
two-form KR field. We can rewrite the action involving 
B as 



S b oc / H A* H - ui L A* u)i 



(A6) 



r i i l 

oc / -dB A* dB + -dB A* u> L + -u L A* dB. 
J 9 3 3 

On variation of this action with respect to B, we have 
the equation of motion, 



d*H = 0. 



(A7) 



Therefore the equation of motion for the KR two-form 
field shows that *H is closed. In other words, at least 
locally, there exists a pseudo-scalar b (the KR axion, or 
sometimes called the universal axion) such that 



H= *db. 



(A8) 



Noting that — *d*d<p = 0<f), we have the equation of mo- 
tion for b, 

Ob = *dH = -- *Tr(R A R) = ^RR (A9) 
6 12 

Varying the action given by Eq. (jAljl with respect to 
the metric we obtain 3 [3] 

- = n 2 a{QH^ x H vXp - g^H^HxpoiAlO) 



-AVa(H Xa ^R 



where G^" is the usual Einstein tensor. Given that the 
equation of motion for the two-form field B allows us to 
write H — *db, we have 



H 



Hup 



We can rewrite the field equation as 



- = n 2 a\2 



rpflU 

1 b 



1 



rXtt(fJ. r> v ) a 



(All) 



(A12) 



where Tt v is the canonical stress-energy tensor for the 
pseudo-scalar field b. We will now show that the last 
term is actually the Cotton- York tensor. 

Using the Bianchi identities for the Riemann tensor, 
we first note that we have the identity, 



(A13) 



With this, it is straightforward to show that 

V CT ([V r 6]e TAQ (^ QA ) 

= 2V CT ([V T 6]/r HrW ) = 2C^, (A14) 

where C^ v is the Cotton- York tensor defined in Eq. ©. 
Choosing a — £ 2 /12 and taking b — * — 9/£ so that in the 
absence of the Cotton- York tensor we regain general rela- 
tivity sourced by a canonical scalar field 9, the equations 
of motion are 



21k' 



C „u = _ K 2 T 



12 



iRR. 



(A15) 
(A16) 



We can see that these field equations are identical to 
Eqs. (j4]) and with vanishing scalar potential. 



3 The sign of the last term in this equation is different in Ref. Q . 
The sign given here makes the divergence of the right-hand side 
vanish as required by the Bianchi identity. 
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APPENDIX B: CALCULATION OF THE 
VECTOR POTENTIAL 

In Lorenz gauge (d^A^ — 0) the Chern-Simons 
Ampere's law, Eq. (|23p . can be written, 



□ 



A- 



1 



-B 



AttGJ, 



(Bl) 



where we have neglected the time variation in 6 in order 
to place m cs inside the d'Alembertian operator. We are 
dealing with a stationary source, and so □ = V 2 . We 
may invert Eq. (|B1[) to obtain 



A 



1 



V x A = G 



J 



rdV. 



(B2) 
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where X is the identity matrix. Multiplying both sides of 



We can write this as 
1 



Vx \A = G 



rdV 



T- (l/m cs )Vx 



the equation by 
A - -^V x V x 4 C ( 7 



we obtain 



-Vx 



J 
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|d J r 
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Noting that V x V x A = — V 2 A in Lorenz gauge, we 
have 



XI 2 A - 



m 2 A = S, 



where 



S = m 2 G ( I - 



-Vx 



J 



rdV. 



(B5) 



(B6) 



We recognize this as the inhomogeneous Helmholtz equa- 
tion. For a rotating homogeneous sphere, the mass cur- 
rent is given by 



J = p[w xr\B(R-r), 



(B7) 



where p is the density, lu is the angular velocity, R is the 
radius, and is the Heaviside step function. 

The most general Green's function for the inhomoge- 
neous Helmholtz equation is 



cos(m cs |f - 7*1) + 7sin(m cs |r - f |) 



Anlr — r 1 



> (B8) 



where 7 is a constant. However, the second term (that is 
proportional to 7) remains constant for (r* — f] <C m~ s , 
implying that the influence of the source does not de- 
crease with distance (for distances r <g; m^. 1 ), which we 
interpret as unphysical. We therefore set 7 = 0. We then 
use multipole expansions for the Green's function, 



cos(m cs |f - 7*1) 



An\f — r'l 

> \ii{^c S r < )yt{m cs r > )Yl m {f')Y lm {f) } 
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l.m 



where je(x) and ye(x) are, respectively, spherical Bessel 
function of the first and second kind, Yi m (r) is a spherical 
harmonic, and the subscript < (>) means the argument 
is the lesser (greater) of r or r'. The solution for A is 
then obtained by integrating, 



.4 



d 3 r'G(r,r')S(r'), 



(BIO) 



where all vectors are expanded in a Cartesian basis. 

The resulting expression for A may be split into a 
general-relativistic and a Chern-Simons term, A = Aqh+ 
Aqs, where 



.4 



GR 



-^i? 3 (rx^)x J* 
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BL 

5r 2 , 



r < R, 

r > R, 
(Bll) 

is the gravitomagnetic vector potential in general relativ- 
ity, and 



AirGpR 3 

TTlraR 



with 

Ci(r) 

C 2 (r) 

C 3 (r) 



r 

m cs r 



[C\(r) lu + C 2 (r) f x lj 
+C 3 {r)f x (f x 3)] , (B12) 



— 2/2 (mo R) ji (m cs r) , 



ml R 2 
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2 R 2 



m cs R y 2 (m cs R)j 1 (m cs r) , 



(B13) 



inside the sphere, and 

. 2i? 3 2R . 

Ci(r) = ——. - H j2(m cs R)yi(m cs r) 

lor a r 

C 2 (r) = m cs Rj 2 (7n cs R)y 1 (m cs r), 
R 3 



C 3 (r) 



5r 3 



7n cs Rj 2 (7n cs R)y 2 (m cs r) 



(B14) 



outside the sphere. We note that this solution for A is 
finite at the origin and continuous across the boundary 
of the sphere, so it produces a finite B at the origin and 
a continuous metric. Taking the curl of this solution for 
A yields the expressions for B given in Section [TV] 

Thus far, we have not discussed any boundary condi- 
tions on the gravitomagnetic field B at the surface of the 
sphere. The field equations for B imply two such bound- 
ary conditions, and we will now prove that the continu- 
ity of A guarantees that these two boundary conditions 
are satisfied. The first boundary condition follows from 
V • B = 0; as in electromagnetism, this condition implies 
that the component of B that is perpendicular to the sur- 
face must be continuous. The second boundary condition 
follows from the Chern-Simons version of Ampere's law: 

VxB - —V 2 B = AnGJ. (B15) 
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Integrating this equation over a surface with vanishing 
area that is perpendicular to the surface of the sphere 
and contains the boundary implies that the components 
of [B + (l/m cs )V x B] that are parallel to the sphere's 
surface must be continuous across the boundary. 

Generally, the continuity of A would not imply conti- 
nuity of its curl. However, our A is a solution to Eq. (|B2jl . 
which may be rewritten as 

A+—B = A GR . (B16) 

m cs 

Since A and Aqr are both continuous across the surface 



of the sphere, this equation implies that B is also continu- 
ous across the surface of the sphere. Furthermore, taking 
the curl of this equation shows that V x B is continuous 
provided that B and V x Aqr are continuous. Taking 
the curl of Eq. (|B1 1|) confirms that V x Aqr is continu- 
ous across the surface of the sphere. Therefore, we have 
shown that the continuity of A implies that both B and 
V x B are also continuous, which guarantees that both 
boundary conditions on B are satisfied by our solution. 
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